Widen et al. Virology Journal 2014, 11:90 
http://www.virologyj.eom/content/1 1 /I /90 



SHORT REPORT Open Access 



PCR detection and analyzis of potentially zoonotic 
Hepatitis E virus in French rats 

Frederik Widen 1 " Florence Ayral 2 , Marc Artois 2 , Ann-Sophie Olofson 1 and Jay Lin 1,3 



□ 



VIROLOGY JOURNAL 



Abstract 

Background: Hepatitis E virus has been detected in a wide range of animals. While Genotypes 1-2 of this virus infect 
only humans, 3-4 can spread from animals to humans and cause sporadic cases of human disease. Pig, and possibly 
also rats, may act as a reservoir for virus. From a public health perspective it is important to clarify the role of rats for 
infection of humans. Rats often live close to humans and are therefore of special interest to public health. Rats live of 
waste and inside the sewage system and may become infected. Reports of hepatitis E virus in rats have been published 
but not from France. The possibility that rats in an urban area in France were Hepatitis E virus infected, with which type 
and relationship to other strains was investigated. This study provides information important to public health and 
better understanding the occurrence of hepatitis E virus in the environment. 

Eighty one rats {Rattus Norvegicus) were captured, euthanized, sampled (liver and faeces) and analyzed by real-time 
RT-PCR's, one specific for Hepatitis E virus in rats and one specific for genotype 1-4 that that is known to infect 
humans. Positive samples were analyzed by a nested broad spectrum RT-PCR, sequenced and compared with 
sequences in Genbank. 

Findings: Twelve liver and 1 1 faeces samples out of 81 liver and 81 faeces samples from 81 captured rats were 
positive in the PCR specific for Hepatitis E virus in rats and none in the PCR specific for genotype 1-4. Comparison 
by nucleotide BLAST showed a maximum of 87% similarity to Hepatitis E virus previously detected in rats and 
significantly less to genotype 1-4. 

Conclusions: This is the first study demonstrating that rats in France carries hepatitis E virus and provide information 
regarding its relation to other virus strains previously detected in rats and other host animals world-wide. Genotype 1-4 
was not detected. 
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Background 

It has been estimated that wildlife is responsible for 72% 
of emergent infectious diseases in humans [1,2]. In order 
to limit the spill over of zoonotic agents, improved wildlife 
pathogen surveillance is required. It is also necessary to 
gain better insight into viral factors for species specificity. 
This is true also for Hepatitis E virus (HEV) a pathogen 
that is present in wildlife, domestic animals and humans 
[3,4]. In developing countries, it causes large scale disease 
outbreaks in humans as well as endemic infections related 
to poor sanitary conditions. In countries with good 
sanitary standards the disease in humans, Hepatitis E, 
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occurs sporadically. Such infections are either acquired 
during travelling to endemic areas or from a domestic 
source in the infected persons home country [3,4]. Wild- 
boars or pigs may constitute such a reservoir [3,4] but 
other species can also play a role. HEV infecting humans 
only or humans and other mammals are taxonomically 
divided in four genotypes [5]. While the large scale epi- 
demics and endemic infections in third world countries 
as well as cases imported from these countries are caused 
by hepatitis E of genotype 1 or 2, infecting humans only, 
the domestically acquired infections in industrialized 
countries are caused by genotype 3 in Europe and North 
America or 4 in East Asia. It is widely believed that pigs 
and wild boars constitute a reservoir for genotype 3 or 4 
for human infections and humans may acquire these 
infections through consumption of undercooked pig, wild 
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boar or deer meat In France the prevalence of anti-HEV 
antibodies in the general population was determined to be 
3.2%, a figure similar to other industrializes countries [6]. 
A survey performed on French pig farms demonstrated a 
65% prevalence of anti-HEV antibodies [7]. Another sur- 
vey on wild boar demonstrated a prevalence of anti-HEV 
antibodies between 7.2 and 22.7% in different geographical 
regions [8]. No figures for Lyon were available but the 
highest value came from Aveyron that is relatively close to 
Lyon. The transmission of HEV between different animal 
species has not yet been clarified. Furthermore it is not 
known what determines the pathogenicity and the host 
range. The current taxonomy is based on a study by Lu 
et al. [5] but has been challenged and additional genotypes 
have been proposed [9]. HEV has been detected in several 
species like for example pig, wild boar, deer, moose, rabbit, 
ferret, mink, rat, poultry and cutthroat trout [3,4,10-15]. 
HEV from pig and wild boar belong to genotype 3 or 4 
and are related. HEV from deer belong to genotype 3 [4] 
and HEV from rabbits is closely related to genotype 3 [9] 
while HEV from moose does not belong to genotype 1-4 
[10]. Avian HEV share only 50% nucleotide identity with 
genotype 1-4 and cutthroat trout HEV is very distantly 
related to the others [14,15]. Rats can be infected with 
seemingly rat specific HEV or with genotype 3 [13,16]. 
The zoonotic potential of rat specific HEV is controversial 
[17,18]. Serological investigations have suggested that 
humans may be infected by rat specific HEV [19]. It is 
likely that genotype 1-4 infecting rats can also infect 
humans. Because R. norvegicus is a synanthropic species, 
humans and rats live in close proximity. Rats are well 
known to transmit pathogens to humans and other 
animals. Their high prevalence and propensity to carry 
pathogens make them a potential reservoir for human 
pathogens including HEV. Indeed, according to previous 
surveys, antibodies against HEV are highly prevalent 
in rats [17] and HEV RNA has been detected in rats 
(R. norvegicus and other rat species) in Germany, USA, 
Vietnam, Denmark, China and Indonesia [13,16,20-23]. 
Differentiation of rat specific HEV from HEV genotype 
3 by serology has been published [19] but is not generally 
available. Most HEV strains found in rats were of the 
rat specific type but genotype 3 RNA has also been 
detected in rats [16]. It is therefore important to further 
clarify the role rats may have as a reservoir for human 
HEV infections. The objectives of this study were to 
investigate if HEV could be detected in French rats in 
an urban environment, to determine the type HEV that 
was present in order to better judge the risk for zoonotic 
infections and finally to determine the relationship to 
other previously detected HEV strains order to create a 
more complete picture of HEV infecting rats. 

The study was part of a survey of infectious agents in 
wild rats (Rattus Norvegicus) conducted in the city of 



Lyon in Center-East France as part of the EU-funded 
Wildtech project. As a part of this survey, liver and faeces 
samples from rats were collected and analyzed for presence 
and characterization of HEV strains. 

Results and discussion 

TaqMan® assay specific for HEV from rat 

Eighty-one rats of the species Rattus Norvegicus were 
caught. No other rat species were observed. Twelve of 
eighty-one (15%) liver samples were positive in real-time 
RT-PCR amplifying HEV from rat [17]. The rats were 
labeled "Rat HEV Ly Id number 2012" and the rats with 
PCR positive liver samples had Id numbers 839, 848, 
867, 873, 874, 877, 879, 880, 882, 883, 888, 894. Testing 
of eighty-one faeces samples resulted in eleven positive 
samples. The positive faeces samples came from rats that 
also had positive liver samples. All positive samples came 
from rats caught in the low income area. Nine (75%) adult 
male and 3 (25%) adult female rats were positive in this 
PCR. No positive juvenile rats were found. 

The sex ratio of the positive rats can be explained by 
the fact that 74% of the rats collected in this area were 
males. Similarly the age ratio of the positive rats (100% 
adults) can be explained by the clear dominance of 
adults (88%) among the collected rats. 

Prevalence and confidence interval of HEV specific for rats 

Fifteen percent (p = 12/81) of the rat tested for HEV had 
a PCR-positive result. Then the 95% confidence interval 
for the prevalence was IC 95 o /o = [7%-22%]. 

TaqMan® assay specific for genotype 1-4 of HEV 

All liver and faecal samples tested by a real-time PCR 
specific for genotype 1-4 [24] were found to be negative 
in this assay. 

Nested PCR of real-time RT-PCR positive samples 

All samples positive by real-time RT-PCR for HEV from rat 
were positive when amplified with this nested PCR [25]. 

Sequence analysis 

Assembly of sequences (DNASTAR Lasergene 8) and 
subsequent BLAST analysis confirmed that the sequences 
represented HEV and were similar to other HEV sequences 
recovered from rats. 

The sequences from the twelve rats displayed highly 
similar sequences. Only one nucleotide position displayed 
synonymous substitution. The HEV strain from rats 874, 
883 and 894 had a "C" at position 4192, as compared to 
the HEV sequence from rat, rat/Mu/0685/DEU2010, 
accession number JN167537.1. Rat 877 had a wobble base 
that varied between "C" and "T" (R) at this position while 
the other sequences had a "T" at position 4192 in the 
cDNA sequence. The multiple sequences deduced from 
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each individual rat sample confirmed this single silent 
mutation. The high level of similarity between these 
sequences can be explained by the fact that all HEV 
positive rats were captured at the same location. Nucleic 
acid BLAST analyses against HEV sequences at NCBI 
gave a maximum identity of 87%, for the entire (100%) 
amplicon sequence, as compared to the German HEV 
strain rat/Mu09/0685/DEU/2010 followed by identities 
between 80% and 86% for other strains from German, 
American, Chinese and Vietnamese rats. The strain 
detected in Vietnam (HEV strain Vietnam- 105, accession 
number JX120573.1), had 82% identity. Two strains from 
USA displayed an identity of 84% (HEV rat/USA/2003, 
accession number JF5 16246.1) and 83% (Hepatitis E virus 
isolate MVZ201020, accession number JQ898482.1). Sev- 
eral Chinese strains displayed identities of 83%. Although 
the difference to other strains from rat, by nucleotide 
sequence BLAST was 12%-20%, the difference by protein 
BLAST was very small (Maximum 98% identity, HEV strain 
KS12/1305, Denmark, accession number AGH06684). The 
high similarity between rat HEV strains is also evident in 
the amino acid alignment (Figure 1). All strains from 
rats collected in this study are 100% identical in the 
amino acid alignment and they differ by a few amino acids 
only when compared to other rat HEV strains. However 
when compared to genotype 1-4 the difference is signifi- 
cantly larger and largest when compared with avian HEV. 

The phylogenetic tree (Figure 2) confirmed that the 
HEV strains detected in this study are most related to 
other Rat specific HEV strains but cluster separately 
from other rat HEV strains. They are most related to Rat 
HEV strains detected in Germany (Figure 2) while the 
strain detected in Vietnam is more distant. 

HEV of genotype 3 has previously been detected in rats 
[14]. These samples were also tested for genotype 1-4 with 
an assay that is regarded as highly sensitive. Because the 
rats were negative for genotype 3 and positive for rat 
related HEV only, they do not pose an evident risk for 
zoonotic HEV infections of humans. The HEV detected 
in these rats is related to other HEV strains from rats. The 
difference compared to HEV genotype 3 is considerably 
larger at the protein level (maximum 74% identity) while 
the similarity is higher when compared with HEV from 
ferrets (maximum 86% identity). The phylogenetic com- 
parison (Figure 2) confirms that detected HEV sequences 
clusters together with other strains from rats and is not so 
distant from the HEV strain from ferret [10]. However, the 
detected strains do not cluster closely with genotype 1-4. 
On the contrary, genotype 1-4 is clustering on a branch 
that is separated from the branch containing rat and ferret 
HEV strains. While it has been clearly demonstrated that 
genotype 1-4 certainly poses a risk for human infections 
there is also a group of mammalian HEV strains that do 
not belong to genotype 1-4. The number of strains in 



this group has been expanding rapidly as an increasing 
number of species have been screened for HEV. While 
some strains like Avian HEV and HEV in cutthroat are 
only distantly related to HEV infecting mammals, there 
are a number of HEV strains that are more similar to 
genotype 1-4 but still sufficiently different not to be 
included in any of these genotypes, like for example rat 
and ferret HEV. These strains are currently not regarded 
as zoonotic. However, further studies of determinants 
of strain specificity, pathogenicity and the relationship 
between different mammalian HEV strains may shed more 
light on light on the zoonotic capability of HEV. 

Methods 

Collection of samples from rats (Rattus Norvegicus) 

Rats provided for this study were trapped for the purpose 
of pest control (trapping agreement n°691810). They were 
captured, euthanized and sampled following the relevant 
ethical and safety rules (animal research agreement n° 
69-020931). The procedure was supervised by the ethical 
committee of the VetAgro Sup and European regulation 
(EU Directive 86/609). 

The necessary sampling size (n«90) was calculated 
using 95% confidence level, a relative accuracy of 50% 
and 10%-30% expected prevalence. 

The survey was conducted in selected areas of Lyon 
during a six month period from October 2011 to March 
2012. The trapping sites were chosen for their abundance 
of rats, as reported by the Hygiene Service of the city and 
their environmental differences relevant for the risk inves- 
tigations. A peri-urban area, a low income area, a public 
garden, a waste treatment plant (WTP) and a waste water 
treatment plant (WWTP), were the five trapping areas. 

The low income area comprised of 619 apartments 
which were distributed on 0.6 km 2 and the traps were 
set in 20% (n = 127) of the dwellings. The public garden 
was situated in the city centre and featured a large pond 
and captive wildlife. In the industrial area, the trapping 
was performed in a WTP and a WWTP. The trapping 
success was 42, 8, 23, 7 and 1 captured Norway rats in 
the low income area, the public garden, the WTP, the 
peri-urban area and the WWTP plant, respectively. A 
total of 81 (56 male and 25 female) rats were captured 
and screened by HEV PCR. Of these 54 were adults (as 
defined by a weight over 100 gram and having sexually 
mature organs) and 27 juvenile. The captured rats were 
labelled "Rat HEV Ly Id number 2012". 

Trapping and sampling of rats 

All rats were captured in small (28 cm x 9 cm x 9 cm) or 
large (50 cm x 15 cm x 15 cm) single catch rat traps. 
The traps were placed on the paths of rats as described 
by people familiar with the area (visitors, employees, 
inhabitants). The distance between each trap depended 
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Figure 1 Multiple sequence alignment of deduced amino acid sequences from nucleotide sequences used for phylogeny (Figure 2). 

The multiple alignment was constructed using MEGA5. It corresponds to the same fragment used for phylogeny (Figure 2). 



on the area to be covered and varied from 1 to 5 m. 
Traps were baited with peanut butter and set at locations 
around the clock for the trapping period. The rats were 
collected each morning. Captured rats were transported 
in the traps to the laboratory, placed inside fume hoods 
and immediately anesthetized using isofluran and eutha- 
nized by cervical dislocation. The species, weight, sex 
and approximate age (juvenile vs adults) of each rat 
was determined by ocular inspection of morphology 
and weighing. Eighty-one rats were aseptically dissected 
and two liver samples of 30 mg each and two faeces 



samples were collected (from the rectum). The samples 
were immediately stored at -80°C prior to shipping to 
SVA in Sweden. 

Extraction of RNA 

Eigthy-one livers samples and 81 faeces samples from 81 
rats were analyzed. The liver samples were cut and the 
fresh surfaces were sampled using a cotton swab which 
was subsequently soaked in 850 ul TE-buffer. The faeces 
samples were diluted 1:4 in TE-buffer. RNA was extracted 
from 90 ul liver sample suspension using a magnatrix 
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Figure 2 Phylogenetic tree depicting the relation between the HEV sequences from French rats to a selection of other HEV sequences. 

Phylogenetic tree of sequences corresponding to a 254 nt long fragment from the nested PCR product of the amplified RdRp fragment. The 
sequence corresponds to nucleotide position 4127 to 4371 of the rat HEV strain rat/Mu/0685/DEU2010, accession number JN 1 67537.1 . The tree 
was constructed by the Neighbor joining method using MEGA 5.05. The tree is depicting the relationship of the French HEV strains from 12 rats, 
here called "Rat HEV Ly, sample number, and Fr 2012" and described in this article, with selected HEV sequences from rat, ferret, avian HEV and 
genotype 1-4. The bar indicates the evolutionary distance as number of base substitutions per site. The bootstrap consensus was generated using 
1000 replicates. Branches corresponding to partitions reproduced in less than 50% bootstrap replicates were collapsed. The percentage of replicate 
trees in which the associated taxa clustered together in the bootstrap test is shown. The tree is drawn to scale. The evolutionary distances 
were computed using the Tamura-Nei method (number of base substitutions per site). The rate variation among sites was modeled with a 
gamma distribution (shape parameter = 6). All positions containing gaps and missing data were eliminated. 



(Magnetic Biosolutions Sweden) or from 50 \i\ faeces 
samples suspension using a MagMax (King Fisher) Robot 

Reverse transcription of RNA and PCR amplification of 
HEV RNA in rat samples 

Eighty-one liver and faeces samples (on from each from 
each rat) were tested by PCR as specified below. 

TaqMan® assay specific for HEV from rat 

Extracted RNA, including positive control R68, was ana- 
lyzed by Taq Man real-time RT-PCR specific for HEV 
from rat, on a Corbett Rotorgene 3000 instrument, as 
previously described [17] but modified according to the 
recommended cycling parameters given in the instructions 
for the Agpath rRT PCR kit (Life technologies, USA) 
and with modified reverse primer (rHEV-R2) and probe 
(rHEV-P2) sequence as recommended by Dr. R. Johne 



(personal communication) to avoid false negative results. 
The PCR target the region 5214-5286 in the rat/Mu/0685/ 
DEU2010 sequence. Briefly; the forward primer rHEV-F 
(5 -TACCCGATGCCGGGCAGT-3 '), the reverse primer 
rHEV-R2 (5 -ATCYACATCWGGGACAGG-3 ') and the 
probe rHEV-P2 (5 -AATGACAGCACAGGCACCGGC 
GCC-3) labelled with 6-FAM at the 5' end and Black 
hole quencher (BHQ) at the 3 ' end was used. The total 
reaction volume was 25 (il per tube including 4.6 \i\ of 
added template. After reverse transcription for 10 min. at 
50°C and inactivation for 15 min. at 95°C, 55 cycles consist- 
ing of 15 sec. at 95°C and 60 sec. at 60°C was performed. 
Fluorescence was collected during the annealing step. 

TaqMan® assay specific for genotype 1-4 of HEV 

Subsequently, all liver and faecal samples, including 
positive control Swe 8, were analyzed by a PCR assays 
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specific for genotype 1-4 as previously described 
[24] with modifications (see below). The PCR target 
approximately position 5060 to 5126 in the rat/Mu/ 
0685/DEU2010 and position 5293 to 5362 in the 
SwX07-El sequence. Briefly, the forward primer JHEVF 
(5'-GGTGGTTTCTGGGGTGAC-3), the reverse primer 
JHEVR (5 -AGGGGTTGGTTGGATGAA-3 ) and the 
probe JHEVP (5 -CCGACAGAATTGATTTCGTCGG 
C-3') labelled with Cy5 at the 5' end and Black hole 
quencher 2 (BHQ2) at the 3' end was used. The total 
reaction volume was 12.5 \A per tube including 3 \A 
added template. The Agpath rRT PCR kit was used for 
the RT-PCR and the RT-PCR program was identical 
with the RT-PCR program specific for HEV from rat. 

Nested PCR of real-time RT-PCR positive samples 

The twelve samples positive by real-time RT-PCR specific 
for rat HEV were amplified with a previously described 
nested PCR targeting a fragment of the RNA dependent 
RNA polymerase (RdRp) in ORF1 [25], Briefly, extracted 
RNA was used as template for a one-step RT-PCR (Qiagen, 
Germany). Primers HEV-cs (5 -TCGCGCATCACMTTYT 
TCCARAA-3 ) and HEV-cas (5 -GCCATGTTCCAGACD 
GTRTTCCA-3 ) were used for the RT-PCR while the 
modified primer HEV-csn mod (5 -TGTTGCCCTGTT 
TGGCCCCTGGTTTAG-3 ') and the primer HEV-casn 
(5 -CCAGGCTCACCRGARTGYTTCTTCCA-3 ') were 
used for the nested step. For the nested PCR, the Platinum 
Taq (Life Technologies) was used. For the one-step 
RT-PCR the reactions were held at 50°C for 30 min., 
and at 95°C for 15 min. followed by 40 cycles of 94°C 
for 30 sec, 50°C for 30 sec. and 72°C for 45 sec. 
followed by final extension at 72°C for 10 min.. For 
the nested PCR the reactions were activated at 94°C 
for 30 sec. followed by 35 cycles of 94°C for 30 sec, 
50°C for 30 sec. and 72°C for 45 sec. followed by final 
extension at 72°C for 10 min. 

Estimation of prevalence 

The estimated prevalence ( p ) in the rat population 
within the sampled area here called "the low income 
area", in Lyon, was the proportion of positive samples. 
The approximate 95% confidence interval (IC 95 o /o ) of the 

true prevalence was given by IC 95 % = p ± 1.96 * with 

"p" the proportion of PCR-positive rats, "q" the proportion 
of PCR-negative rats and "n" the sample size. 

Sequencing of HEV positive samples 

The 12 samples successfully amplified by nested PCR 
were sequenced from different PCR products on average 
4-5 times per product using the inner primers from the 
nested PCR. The sequencing was performed at an 
Applied Biosystem 3130x1. 



Sequence analysis 

Multiple alignments of sequences were performed using 
DNA-STAR Lasergene 8 and a 254 nt fragment was 
selected for further analysis. This fragment was chosen 
because it was the largest that had good sequence 
coverage for all 12 rat samples. Consensus sequences 
for individual rat samples were created with DNASTAR 
Lasergene 8. These consensus sequences were analyzed 
using the nucleotide BLAST program at NCBI as well 
as by Mega 5.05. 

Phylogenetic and molecular evolutionary analyses were 
conducted using MEGA version 5.05. The phylogenetic 
tree was deduced using Neighbour joining with Tamura 
3 parameter and gamma parameter 0.3. Bootstrap analysis 
was performed with 1000 repeats. 
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